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ABSTRACT 

The results of a study of the AGB phase of stellar evolution are presented. 
Abundances have been determined for Fe, C, O, the light s-process elements, Y 
and Zr, the heavy s-process elements, La and Nd, and the r-process element, 
Eu. The expected relationship between enhanced C, increasing C/0 ratio and 
enhanced s-process elements has been quantified. Results are presented to provide 
observational data with which to compare theoretical predictions. 

The results in this paper confirm previously suggested relationships between 
C, C/0 and s-process element enhancements. It is seen that AGB stars show 
G/0 ratios from G/O~0.4 to 1.0, while G enhancements lie between [G/Fe]=0.1 
to 0.9 dex. Enhancements of s-process elements are as much as [s/Fe]~1.0 dex 
for the stars in which G is also greatly enhanced. 



Subject headings: abundances - nucleosynthesis - stars: AGB and post-AGB 



- 2 - 



1. Introduction 

The analysis of the surface abundances of asymptotic giant branch (AGB) stars offer 
tests for the operation of the third dredge-up (TDU). 

The surface of AGB stars can be enriched in carbon as a direct result of the TDU 
phenomenon (Lattanzio and Karakas, 2001). Stars are formed with a C/0 ratio <1 but at 
the end of the AGB the stars can show a C/0 ratio >1 in the envelope as a result of this 
mixing process that pollutes the envelope with carbon-rich material. The stellar mass plays 
a critical role in determining the enhancement of both carbon and the s-process elements, as 
it directly affects both the source of neutrons in the stellar interior and the efficiency of the 
TDU episodes (see Busso et al., 1999, and references therein). The final C/0 ratio is also 
affected by the mass loss history although a discussion of this is not given here. Thorough 
details of observations of AGB stars are vital to enable comparison with the enrichments 
predicted by numerical stellar nucleosynthetic modeling. Stars can either enrich themselves 
via their own internal TDU, in which case they arc known as intrinsic carbon stars, and 
display technetium in their spectra, or arc enriched from transfer of carbon-rich material 
from a binary companion, in which case they are known as extrinsic carbon stars and any 
(radioactive) technetium will have decayed into a stable isotope of ruthenium. 

The carbon content of the stellar envelope is expected to increase with changing spectral 
type along the AGB, from M to MS to S to C, where M stars show C/O~0.5 while C stars can 
show C/0>1. This TDU enrichment process also acts to increase the envelope abundances 
of s-process elements. Light s-process (Y and Zr) and heavy s-process (La and Nd) elements 
discussed in this study are hereafter referred to as 'Is' and 'hs' respectively. 

The spectra of M stars are classified by strong absorption bands of TiO and by the large 
number of strong metaUic fines blueward of 4000A(Jaschek and Jaschek, 1987). The strength 
of the TiO bands increases rapidly with advancing subtype, and continuum fitting becomes 
difficult. Near the end of the M classification, bands of VO become strong compared to TiO. 
MS stars are placed between the separate M and S classifications, and display molecular 
bands of both ZrO and TiO, although the latter still dominates. 

S stars have a C/0 ratio approaching unity. They are late- type stars showing distinct 
bands of ZrO in the blue and visual. The easiest way to distinguish S stars is the compar- 
ison between strengths of the ZrO and TiO bands, with ZrO taking over as the S subclass 
increases. 

SC stars are thought to fie between S and C type stars and show characteristics of both 
classes. There is some observational evidence that SC stars can be formed in two different 
ways; either through the continuous dredge-up of C to the envelope during the AGB phase. 
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or due to the hot bottom burning phase in massive AGB stars. In either case, the phase 
during which C/0 ~ 1.0 is very short, which may explain why the number of confirmed SC 
stars is so low (Guandalini and Busso, 2008). They are classified by their C/0 ratio of 
with most of the C and O locked up in CO. They exhibit enhanced heavy metal lines, and 
bands of CN are prominent. They also display very strong NaD lines and weak but clear 
C2 and ZrO bands. Due to the balance of oxygen and carbon, all oxide bands and carbon 
molecules are weak and the spectra are full of atomic lines, making them easier targets for 
abundance analysis. However, only a few in-depth studies on heavy element abundances 
have been undertaken (Catchpole, 1982; Abia and Wallerstein, 1998). 

C stars show very crowded spectra due to strong molecular bands and low temperatures. 
Their C/O ratios are very similar to those shown by SC stars. Again, only a few studies on 
heavy element abundances have been undertaken (Utsumi, 1970; Abia et al., 2001; Abia et 
al, 2002). 

In order to compare observed s-process enhancements with those predicted by theory, it 
is desirable to have a set of results along the AGB phase that are internally self-consistent. 
Relatively few studies have been devoted to large samples of AGB stars, usually studying 
only one or two spectral classes at a time. 

A sample of 22 M and MS stars have been analysed by Smith and Lambert (1985, 1986). 
These studies used high-resolution, infra-red spectra to obtain abundances for C, N, O, Is 
and hs elements. They found s-process enhancements in seven of the eleven M stars but no 
measurable enhancements in four of the eight MS stars. 

The most detailed study of S and SC stars has been undertaken by Abia and Wallerstein 
(1998) using high-resolution spectra over the optical wavelength range. This study of S and 
SC stars derived abundances for C, N, O and a large number of heavy elements. They found 
considerable s-process elements enhancements in all S and SC stars, of values as much as 
[s/Fe] — hl.4 dex. 

Both of these studies measured equivalent widths for atomic lines, including s-process 
lines, which are often weak and are affected heavily by molecular blending. 

This present study of eleven stars has eight stars in common with these previous studies, 
and includes three stars not previously analysed. It provides a self-consistent set of results 
over a range of spectral classes, M, MS, S and C using high-resolution (R~40,000), high 
signal to noise (S:N>100) spectra over an optical wavelength range of 4500 to 75OOA. The self- 
consistency of these results minimises all discrepancies due to differing methods or systematic 
offsets between studies. It allows for an internally accurate abundance analysis of stars 
relative to each other along the AGB phase. Observations and analysis of this sort are of 
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great importance to compare to the enhancements predicted by stellar models in order to 
refine various variable model parameters. Full details of the method of analysis are discussed 
in Section 2. The results of the abundance analysis are presented in Section 3, before a 
discussion of both specific results and general trends in Section 4. A brief conclusion is given 
in Section 5. 



2. Observations and Analysis 

The eleven field AGB stars analysed are listed in Table [H with previously published spec- 
tral types (Smith and Lambert, 1986), recent magnitudes (SIMBAD databas^ and 2MASS 
survejo), and calculated bolometric magnitudes, variability type (Samus and Durlevich, 
2004) and previous detection of the radioactive element, technetium (Smith and Lambert, 
1988). Bolometric magnitudes were calculated using the bolometric correction approxima- 
tion from Alonso et al.(1999) and absolute K magnitude s using SIMBAD's a pparent mag- 



nitudes and distances from the HIPPARCOS catalogue flHIPPARCOSl 119971 ). Magnitudes 
are not given for HD286340 due to its extreme variability preventing reliable measures of its 
brightness. 

Observations were taken with the Im telescope and HERCULES spectrograph at Mt 
John University Observatory, New Zealand, over 2004 September, and the Anglo-Australian 
Telescope and UCLES at Siding Springs over 2004 August. Standard reduction procedures 
in FIGARO were followed for all observations. 

The abundance analysis for these field stars included using atomic line lists from Kuruc^ 
with refined oscillator strength values for the s-process lines via a reverse solar analysis. This 
reverse solar analysis was done with a solar model of Te//=5770K and log g=4.44. A similar 
reverse analysis was also undertaken using Arcturus as this red giant more closely represents 
the atmospheres dealt with in this research and lines too weak for detection in the sun are 
often seen at these cooler temperatures. The Arcturus model adopted was Te//=4300K, log 
g=1.50 and [Fe/H]=-0.50. Due to the fact that at these cooler temperatures and in these 
C enhanced stars, molecular features become very strong, all molecular features of CO, CN, 
CH, ZrO, VO, and TiO were included. With the exception of TiO line data, which was 
taken from Plez (1998), all molecular line data was taken from Kurucz's linelists. This study 



^http: / /simbad.u-strasbg.fr 
^http:/ /irsa. ipac.caltech.edu 
•^http: / /kurucz. harvard.edu 
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does not employ equivalent width measurement for any lines other than iron due to strong 
molecular blending and difficult continuum placement. All abundances are obtained using 
the spectrum synthesis program MOOG (Sneden, 1973). The Fe lines used for equivalent 
width analysis to determine the atmospheric parameters (see Section 2.1) are listed in Table 
|2] while the oxygen and heavy element lines used for spectrum synthesis analysis are listed 
in Table [31 This spectrum synthesis program adopts solar abundances from Anders and 
Grevesse (1989). Of specific note are the carbon, nitrogen and oxygen abundances, which 
were taken as log(e)C=8.56, log(e)N=8.05 and log(e)0=8.93. 



2.1. Atmospheric Parameter Determination 

Due to severe molecular blending and difficult continuum placement in these stars, it 
becomes increasingly difficult to find unblended Fe I and Fe II lines for equivalent width 
measurement with which to determine the atmospheric parameters. However, this technique 
was adopted for all stars in an attempt to obtain a spectroscopic effective temperature and 
gravity. With a total of around 28 Fe I lines and 8 Fe II lines, the usual comparison between 
excitation potential and derived abundance was used to determine the effective temperature 
while the balance between abundances derived from neutral and ionised iron lines was used 
to refine the gravity of the star. The microturbulent velocity was obtained by measuring Fe 
I lines and ensuring that the derived abundances be independent of equivalent width. 

The final spectroscopic effective temperatures ranged from 3650K to 4000K. A pho- 
tometric temperature was also calculated using the colour-index calibration published by 
Alonso et al. (1999), assuming an [Fe/H] of 0.0 for this calibration (see Table H]), which 
ranged from 3640K to 3970K. With three exceptions, HD7351, HD49331 and HD64332, the 
spectroscopic and photometric temperatures agree to within 50K. All gravities fell between 
the expected values of log g=0.5 to 1.5. For one star, IID286340, no reliable B-V values could 
be found due to the star's extreme variability. In this case, typical values of Tejj=3500K, 
log g=1.0 and ^=1.50 km s~^were chosen. All atmospheric parameters derived are shown in 
Table |H Model atmospheres are taken from the MARCS model atmosphere grid (Gustafs- 
son et al, 2008). C-rich MARCS model atmospheres were adopted for those stars in which 
C/O~1.0 (HD 30959, HD 64332 and HD 286340). 
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2.2. 



C and O Line Detail 



In addition to tlie s-process element abundances, tlie C and O abundances are of rel- 
evance along the AGB phase. Given that the C/0 ratio changes along the AGB, due to 
the dredging up of C after the thermal pulses, it is expected that a correlation should exist 
between [C/Fe], C/0 and [s/Fe]. In order to quantify this relationship, C and O abundances 
were determined for these field stars. 

The C abundance was determined from the C2 bands in the 5060 to 5100A region. The 
nature of C2 results in extreme sensitivity to any change in C abundance, making it a reliable 
molecule to use in deriving any C enhancements. The C2 molecular line parameters were 
taken from Kurucz's linelists. 

The O abundance was obtained from spectrum synthesis of the forbidden O lines at 
6300 and 6363A. The oscillator strength of the forbidden line at 6300A is well known, and 
a value of log(g)=-9.750 was adopted (as discussed in Lambert, 1978). The 6363A O line 
is very weak and should not be used for abundance determination on its own. However, in 
this study it was used as a check of the abundance derived from the 6300A line wherever it 
was detectable. The oscillator strength adopted was log(g)=-10.25 (Lambert, 1978). These 
values differ from those supplied in Kurucz's line list by only -0.07 and -0.05 respectively. 



All abundance results for this sample of field AGB stars are shown in Tables [5] and [61 
Specific results of interest are discussed further below. 



Iron abundances were obtained for all stars for which spectroscopic atmospheric param- 
eters could be determined. The solar iron abundance adopted for all abundance analyses 
was log(e)Fe = 7.52, close to the meteoritic abundance published by Anders and Grevesse 
(1989). For HD286340, for which standard parameters were adopted, an iron abundance 
of [Fe/H]=0.0 was assumed. [Fe/H] values ranged from -0.26 dex to +0.17 dex. Thus all 
stars analysed can be considered of solar or near-solar metallicity. All stars are members 
of the disk population, and the observed spread in [Fe/H] in the sample results from both 
uncertainties and genuine differences. Of the eight stars in common with previous studies 
(Smith and Lambert, 1985; Smith and Lambert, 1986; Abia and Wallerstein, 1998), all but 



3. 



Results 



3.1. 



Fe abundance 
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one [Fe/H] value agrees with previous reported values to within 0.10 dex. 

As a check of the non-LTE effects on low-excitation potential neutral lines the iron 
abundance was plotted against effective temperature, as seen in Figure [TJ No real trend is 
seen for [Fe/H] as a function of temperature. The Fe I lines used for effective temperature 
determination have medium to high excitation potentials, x = 2.5 — S.OeV, and are thought 
to be less affected by dramatic non-LTE effects. 

3.2. C and O abundances 

As mentioned earlier, the abundances of C and O in AGB stars are of interest when 
comparing with theoretical predictions. It is expected that the C abundance will increase 
with recurring thermal pulses as fresh carbon is dredged to the surface. Consequently, the 
C/0 ratio changes as the star moves up the AGB, from C/O~0.5 in M stars to C/O~1.0 
in C stars. Figure [3] shows examples of the C2 region (at 5140A) and the O I forbidden 
line (at 6300A), used to derive the C and O abundances. Syntheses with no C or O present 
are shown as the short dashed lines. The other synthetic spectra show [C/Fe] at varying 
enhancements and [O/Fe]=0.00, +0.10, +0.25 and +0.40. 

As can be seen in Tables O and O most stars show enhancements in C from [C/Fe]=+0.11 
to +0.79 depending on spectral type, as expected. Oxygen was also observed to be either 
near-solar or slightly enhanced, from [O/Fe]=0.00 to +0.30, again depending on spectral 
type. Generally speaking, the two elements are enhanced by similar amounts in any one 
star, with stars greatly enhanced in C also showing slight enhancements in O. The C/0 
ratios are displayed in Figure [21 with C/0 shown as a function of the frequency distribution 
among the sample. The C/0 ratio is an indication of the quantity of carbon brought to the 
surface by the third dredge- up process. As can be seen from Tabled the M stars have C/0 
ratios ranging from C/O=+0.47 to +0.79 with a mean of C/O=+0.59. Table [6] show the S 
and SC stars values ranging from C/O=+0.83 to +1.32 with a mean of C/O=+1.00. There 
is quite clearly a distinct spread in the C/0 values in this sample. 

3.3. s-process element abundances 

The M stars show small or no enhancements of s-process elements, Y, Zr, La and Nd. 
Four M stars give mean values of [ls/Fe]=+0.08 and [hs/Fe]=+0.09, while the other outlier, 
HD7351(M3) gives more enhanced values of [ls/Fe]=+0.45 and [hs/Fe]=+0.46. These four 
M stars can be said to show minimal s-process element enhancements. This is consistent 
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with theory, in that there is not expected to be any process in M stars that would cause 
enhanced s-process element abundances as there have not been any thermal pulses this early 
in the AGB phase. The star HD7351 also displays a higher C/0 ratio, which suggests it 
is a more evolved spectral class, discussed further in later sections. Of the eight stars in 
common with previous studies (Smith and Lambert, 1985; Smith and Lambert, 1986; Abia 
and Wallerstein, 1998), most s-process element abundances agree to within 0.18 dex. The 
greater the observed enhancement the more the divergence in agreement between studies, 
although abundances in this study never disagree with previous studies by more than 0.25 
dex. There also appears to be no systematic offsets in this analysis as abundances in this 
study are equally distributed both higher and lower than previous published values. 

Figure m shows examples of a Zr I (at 6140A) and Nd II line (at 516lA), used to derive 
the Zr and Nd abundances. Syntheses with no Zr or Nd present are shown as the short 
dashed lines. The other synthetic spectra show [X/Fe]=0.00 and enhancements of +0.5 and 
+ 1.0 dex. It is clear that the M star (3a, 3b) shows no Zr or Nd enhancements while the C 
star (3c, 3d) shows increasingly strong Zr and Nd lines. Figures [5] and [6] show schematics of 
slightly larger regions and can be used to compare the relative strengths of Zr and Nd lines 
in a K giant (top: HD 27371) compared to M (HD 49331), S (HD 49368) and C (bottom: 
HD 30959)) stars. 

There is a definite separation, marked by an increase in the s-process elements in the S 
stars. The S and C stars show mean values of [ls/Fe]=+0.76 and [hs/Fe]=+0.65. A single 
star, HD216672 (S5) has considerably lower values of [ls/Fe]=+0.39 and [hs/Fe]=+0.39 dex. 
With its lower s-process element enhancements and lower C/0 ratio, it and HD7351 lie 
between the two extreme groups of M and C stars. 

3.4. Uncertainties on derived abundances 

All abundances quoted in this paper carry uncertainties, some of which are more signif- 
icant than others. There are three main sources of uncertainty in these derived abundances, 
discussed here in detail. 

The main source of uncertainty is the spread in derived abundances from multiple lines 
of the same element. While, in theory, this should be small, different lines of the same 
element can give different abundances, by as much as ~0.2 dex. This can be due to incorrect 
log gf values, although in this study a reverse analysis using Arcturus was done for all lines 
used to avoid precisely this problem (see Section [2] for details) . Non-LTE effects, blending 
and data quality may also contribute to this problem. It is this value that is quoted as 
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the a value in Tables [5] and [61 Where abundances are obtained from only one line (C, O, 
Eu and occasionally s-process elements) the a value is the uncertainty from a single line 
measurement (discussed below). 

The second source of uncertainty is the accuracy with which it is possible to determine 
an abundance from any particular individual line. This uncertainty can vary markedly 
depending on the line strength, any blending or crowding in the region and the overall 
quality of the data. In this study it is estimated that this value should be no more than 
~ O.lSdex. 

The final source of uncertainty is the difference in abundances obtained using different 
atmospheric parameters in the stellar models. This depends solely on how well the atmo- 
spheric parameters can be constrained. Table [7] shows the dependence of derived abundances 
on the chosen atmospheric parameters. These values were obtained by choosing equivalent 
width values which gave abundances representative of those obtained via spectrum synthe- 
sis and then running an abundance analysis on these equivalent widths with four different 
models, with ATe//=+250K, Alog g=+0.3 and A^=+0.25km s^^. These A values repre- 
sent higher uncertainties on the atmospheric parameters than the ones estimated in this 
study (see Table H]). It can be seen that all atmospheric parameters influence different lines, 
stressing the importance of ascertaining a correct model atmosphere at the outset. 

Two elements which merit specific discussion are C and O. Given the importance of 
these two elements to AGB modeling, it is necessary to provide a realistic estimate of the 
uncertainty involved in these abundance derivations. The abundances for C and O are 
derived from one region or line only so the uncertainty associated with these abundances 
depends on how well that particular feature can be measured, as opposed to an abundance 
spread over various features. It is estimated that the C and O abundances can be derived 
to within ±0.15 dex on these particular features (examples shown in Figure [3]). This has a 
dramatic effect in the uncertainty in the quoted C/0 ratios. Taking into account the 0.15 
dex uncertainty in both C and O, it results that the uncertainty on the final C/0 ratio is 
±0.30 dex. While this is a large quoted uncertainty it is hoped that, in reality, the C/0 value 
has an uncertainty less than this. Data of higher quality and measurement of more features 
of C and O, including molecular features out of wavelength range of this investigation, would 
refine this value further. In addition, these values are absolute uncertainties and the relative 
uncertainties between the stars should be smaller due to the consistent data sets and analysis. 
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4. Discussion 

4.1. General trends in AGB stars 

Generally, it is expected that as the star moves up the AGB the [C/Fe] value and C/0 
ratio increase due to the C^^ brought to the surface by the third dredge-up process. Another 
important by-product of this phenomenon is the expected surface enhancements of the s- 
process elements. In an attempt to quantify this relationship, this research examines the 
[Is/Fe] and [hs/Fe] values as a function of increasing AGB spectral type (Figure [7]), [C/Fe] 
value (Figure [8]) and C/0 ratio (Figure [9]). 

Figure [7] shows the enhancements of the light s- (as filled circles) and heavy s- process 
elements (as open circles) for each individual star as the AGB phase progresses. The M stars 
are clearly minimally enhanced, and as the spectral type advances the enhancements of the 
light and heavy s-process elements clearly increases. However, this figure relies heavily upon 
spectral classification, which can often be difficult or ambiguous. 

A more logical way of quantifying s-process enhancements is through the [C/Fe] and 
C/0 values. Figure [H] shows the [Is/Fe] (as filled circles) and [hs/Fe] (as open circles) values 
as a function of the derived [C/Fe] value. In general, it is clear that as the [C/Fe] value 
increases s-process enhancements are seen. All stars seem to follow this trend between 
[C/Fe] and [s/Fe]. However, there is some indication that the s-process elements may be 
enhanced first, followed by the carbon enhancement. This is most clearly seen in Figure [HI 
in which there is a steep rise in [Is/Fe] and [hs/Fe] at [C/Fe] = 0.10 to 0.50 and then a plateau 
region from [C/Fe] = 0.50 to 1.0, at which the s-process elements have an approximately 
constant enhancement of ~1.0. Although this generalization is based is on a relatively small 
number of objects, it may provide some observational test data for the modelers of the TDU. 
It should be noted that this sample most likely includes a mix of stellar masses and this is 
a critical parameter for the efficiency of the TDU episodes and, consequently, the C and 
O abundances. It is very difficult to obtain accurate masses for field stars, and thus no 
estimate of stellar mass has been attempted here. This problem can be avoided by analysing 
AGB stars in a globular cluster, where masses are more accurately determined although few 
studies have been published (Wylie et al., 2006; Yong et al., 2008). 

Figure [9] shows the [Is/Fe] (as filled circles) and [hs/Fe] (as open circles) values as a 
function of C/0 ratio. This ratio is also an indication of the amount of dredge-up occurring 
in the star and is expected to increase from ~0.5 at the beginning of the AGB phase to ~1.0 
at the end of the thermally pulsing AGB stage. Results from several stars in previous studies 
(Smith and Lambert, 1985; Abia et al. 2002) that are not in common with this paper are 
also included as squares and triangles respectively. 
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As can be seen in Figure [H generally all stars follow the expected trend of increasing 
[s/Fe] with increasing C/0 and the spread does cover the expected range of C/O=0.5 for 
M stars to C/O>1.0 for C stars. The plateau effect at [Is/Fe] and [hs/Fe] = +1.0 is also 
evident. Figure also shows how well this current study agrees with trends seen in previous 
studies. The stars from Smith and Lambert (1985), shown as squares, and the stars from 
Abia et al. (2002), shown as triangles, seems to show no marked discrepancies from the 
trend observed in this study. The inclusion of these 17 stars, which results in a larger overall 
sample, strongly supports the general conclusions drawn from the smaller sample studied in 
this paper. 

Another result that Figure [9] suggests is that as the C/0 ratio increases the difference 
between the [Is/Fe] and [hs/Fe] enhancements gets larger. This can be confirmed by looking 
at the individual [hs/ls] values as a function of C/0. For the three stars with C/0 > 0.95 
their average [hs/ls] value is -0.22, whereas the average [hs/ls] for all other stars is -0.02. 
This may suggest that initially, while the s-process and carbon enhancements are minimal 
they are all enhanced by the same amounts. However, as the enhancements increase in value, 
the two s-process element peaks separate and begin to be enhanced at different rates, with 
the hght s-process peak being favoured. For near-solar metallicity it is predicted that the 
light s-process elements are more enhanced than the heavy s-process elements (see Busso 
et al. (2001) and references therein), and earlier studies have also suggested this trend 
observationally (Smith and Lambert, 1985; Abia et al. 2002). This prediction is now further 
confirmed by these present results. 



4.2. r-process, Eu 

As can be seen in Tables [5] and El three stars, HD44478 (M3), HD112300 (M3) and 
HD216672 (S5) show slight deficiencies in the r-process element, Eu, while three others, 
HD49368 (S3), HD64332 (S4) and HD286340 (SC7) show possible small enhancements in 
Eu (of only [Eu/Fe]~0.1dex), which are not predicted by the theoretical modelling. These 
enhancements are small and within uncertainties the [Eu/Fe] ratio is compatible with solar 
ratios. However, if these enhancements are genuine, and not a product of uncertainty on 
measurements, this raises an interesting issue. The AGB evolutionary phase is not expected 
to produce r-process elements in any way, so it is not expected that AGB stars would show 
Eu enhancements due to their stellar evolution. These observed enhancements may be more 
likely to come from the primordial enrichment of the material from which the stars formed 
or external pollution from other sources. 

If this Eu enhancement is genuine, it is possible that these three stars have undergone 
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pollution from material rich in r-process elements. This may have involved extrinsic pollution 
from supernova remnants or other ejecta in which r-process element formation has previously 
occurred. This may explain any observed Eu enhancement in these stars. It may also help 
to explain the observed Nd enhancement in these stars, which in each of these three stars 
shows a greater enhancement than the other heavy s-process element, La. Nd has five stable 
isotopes, four of which are formed by both the s- and r-process and one of which is r-process 
only. If these two stars had undergone pollution from r-process element rich material, then 
it is possible that the source of this pollution is responsible for the observed enhancement of 
Nd. La, which has one stable isotope, can also be produced by both the s- and r-process, thus 
providing more support for this theory. The large enhancements of La, Nd and the slight 
enhancement in Eu observed in these stars may not be due solely to AGB nucleosynthesis, 
but partially due to pollution from previously r-processed matter. 

5. Conclusion 

In this paper the results of an internally accurate and self-consistent abundance analysis 
of AGB stars are presented. The relationship between [C/Fe], C/0 and s-process element 
enhancements is quantified. All general trends seen in this small sample of stars are strongly 
supported by those found in previous studies. While absolute uncertainties on the crucial 
elements C and O are reasonably large at ±0.30 dex, the relative uncertainties should be 
smaller due to consistent data sets and analysis methods, thus supporting the strength of 
the results from this internally consistent analysis. 
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Table 1. Basic stellar quantities. 



Star 


Spectral Type^ 


B 


B 


K 


MBol 


Variable Type'' 


Tc detection'^ 


HD7351 


M3 / S 


8.07 


6.38 


1.64 


-3.02 


SR 


no 


HD30959 


M3 / S 


6.53 


4.75 


-0.66 


-3.87 


SRB (30 days) 


yes 


HD35155 


S 


8.64 


6.87 


2.14 


-4.45 


EA/GS/WD 


no 


HD44478 


M3 


4.53 


2.91 


-1.86 


-3.31 


LB 




HD49331 


M2 / S 


6.91 


5.11 


0.56 


-4.12 


susp 


no 


HD49368 


S3 


9.53 


7.78 


2.46 


-3.57 


SRB 


no 


HD64332 


S4 


9.34 


7.64 


2.31 


-2.42 


LB 


yes 


HD102212 


Ml 


5.58 


4.05 


0.16 


-2.10 


SRB 




HD 11 2300 


M3 


4.96 


3.38 


-1.19 


-2.37 


susp 




HD216672 


S5 


8.21 


6.47 


1.04 


-3.42 


SRB (50 days) 


yes 


HD286340 


SC7 










SRB (370 days) 





HSmith and LambertI lll986h 

ISamus and DurlevichI ||2004| V SR: semi-regular, SRB: semi-regular with poorly defined period, 
EA/GS/WD: eclipsing binary with both giant and white dwarf component, LB: slow irregular, susp: 
only suspected variable, type not confirmed. 

iSmith and LambertI lll988h 
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Table 2. Line list of all Fe I and Fe II lines used to derive abundances via equivalent 

width analysis. 



Species 


Wavelength 


X 


log gf 


Species 


Wavelength 


X 


log gf 




(A) 


( eV) 






(A) 


(eV) 




Fe I 


5044.21 


2.85 


-2.150 




5859.56 


4.55 


-1.080 




5068.77 


2.94 


-1.230 




5862.33 


4.55 


-0.600 




5074.75 


4.22 


-0.200 




5929.68 


4.55 


-1.410 




5159.06 


4.28 


-0.820 




6065.48 


2.61 


-1.530 




5215.18 


3.27 


-0.870 




6094.12 


4.65 


-1.640 




5618.63 


4.21 


-1.260 




6096.66 


3.98 


-1.830 




5624.02 


4.39 


-1.330 




6107.91 


4.14 


-0.191 




5633.95 


4.99 


-0.120 




6226.74 


3.88 


-2.220 




5635.82 


4.26 


-1.740 




6240.64 


2.22 


-3.388 




5662.50 


4.12 


-0.570 




6380.74 


4.19 


-1.401 




5686.51 


4.55 


-0.450 


Fe II 


5197.56 


3.23 


-2.100 




5691.50 


4.30 


-1.370 




5234.62 


3.22 


-2.050 




5705.47 


4.30 


-1.360 




5534.85 


3.24 


-2.740 




5717.83 


4.28 


-0.980 




6084.10 


3.20 


-3.908 




5731.76 


4.26 


-1.150 




6238.39 


3.89 


-2.630 




5806.73 


4.61 


-0.900 




6247.55 


3.89 


-2.516 




5809.22 


3.88 


-1.690 




6432.68 


2.89 


-3.808 




5852.22 


4.55 


-1.180 




6456.39 


3.90 


-2.275 
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Table 3. Line list of all other lines used to derive abundances. 



Species 


Wavelength 


X 


log gf 


Species 


Wavelength 


X 


log gf 




(A) 


(eV) 






(A) 


(eV) 




O I 


6300.30 


0.00 


-9.75 




6140.46 


0.52 


-1.41 




6363.78 


0.02 


-10.25 




6143.18 


0.07 


-1.10 


Y I 


5526.72 


2.00 


-0.65 


Zr II 


5112.28 


1.66 


-0.59 




6435.05 


0.07 


-1.02 




5350.09 


1.83 


-0.94 


Y II 


5087.43 


1.08 


-0.17 




5350.36 


1.77 


-1.18 




5119.12 


0.99 


-1.36 


La II 


5114.56 


0.24 


-1.06 




5473.39 


1.74 


-1.02 




5122.99 


0.32 


-0.93 




5544.62 


1.90 


-1.09 




5797.60 


0.24 


-1.51 




5546.03 


1.75 


-1.10 




5805.77 


0.13 


-1.61 


Zr I 


4805.89 


0.69 


-0.57 




5808.31 


0.00 


-1.56 




4815.06 


0.65 


-0.53 


Nd II 


5161.71 


0.74 


-0.98 




4815.64 


0.60 


-0.13 




5165.13 


0.68 


-0.06 




4828.06 


0.62 


-0.64 




5167.92 


0.56 


-0.98 




4809.48 


1.58 


-1-0.16 




5795.15 


1.26 


-1.13 




6127.48 


0.15 


-1.06 




5804.00 


0.74 


-0.53 




6134.57 


0.00 


-1.28 


Eu II 


6645.13 


1.38 


4-0.40 
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Table 4. Atmospheric Parameters for Field Stars 



Star 


Ph. T,^^ 
(K) 


Sp. T^ff 
(± 200K) 


log g 
± 0.2 


{± 0.25 km s-l) 


HD7351 


3740 


3650 


0.5 


1.50 


HD30959 


3650 


3600 


1.0 


1.40 


HD35155 


3670 


3700 


1.0 


1.30 


HD44478 


3850 


3900 


1.0 


1.80 


HD49331 


3640 


3500 


0.5 


1.25 


HD49368 


3700 


3700 


1.0 


2.00 


HD64332 


3750 


3600 


1.0 


1.50 


HD102212 


3970 


4000 


1.5 


1.40 


HD112300 


3890 


3900 


0.5 


1.10 


HD216672 


3700 


3700 


1.0 


1.20 


HD286340 




3500 


1.0 


1.50 
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Table 5. Element Abundances for field M and MS stars 



X/H 


HD7351 


a 


HD44478 




HD49331 


c 


HD102212 




HD112300 


a 


Fc I 


0.05 


0.08 


0.00 


0.09 


0.15 


0.29 


-0.04 


0.12 


0.02 


0.05 


Fe II 


0.02 


0.10 


0.00 


0.08 


0.22 


0.06 


-0.04 


0.08 


0.00 


0.06 


[X/Fe] 






















C 


0.28 


0.15 


0.15 


0.15 


0.13 


0.15 


0.18 


0.15 


0.11 


0.15 


o 


0.01 


0.15 


0.00 


0.15 


0.03 


0.15 


0.08 


0.15 


0.07 


0.15 


Y I 


0.29 


0.15 


0.04 


0.03 


0.16 


0.03 


0.14 


0.15 


-0.01 


0.15 


Y II 


0.39 


0.13 


0.06 


0.02 


0.10 


0.03 


0.16 


0.10 


-0.01 


0.10 


Zr I 


0.50 


0.24 


0.06 


0.03 


0.16 


0.06 


0.26 


0.09 


-0.08 


0.02 


Zr II 


0.51 


0.21 


0.03 


0.07 


0.11 


0.07 


0.29 


0.10 


-0.01 


0.07 


La II 


0.43 


0.07 


0.05 


0.04 


0.08 


0.01 


0.27 


0.04 


-0.02 


0.15 


Nd II 


0.49 


0.07 


0.09 


0.03 


0.07 


0.15 


0.19 


0.07 


-0.01 


0.09 


Eu 11 


CO,') 


0.10 


-0.11 


0.10 


-0.01 


0.10 


0.02 


0.10 


-O.O:-! 


0.10 


C/O 


0.79 




0.60 




0.54 




0.54 




0.47 




Is 


0.45 




0.04 




0.08 




0.22 




-0.02 




hs 


0.46 




0.07 




0.07 




0.23 




-0.02 




hs/Is 


0.01 




0.02 




-0.01 




0.01 




0.00 
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Table 6. Element Abundances for field S and SC stars 



X/H 


HD30959 


c 


HD35155 


c 


HD49368 




HD64332 


c 


HD2 16672 




HD286340 


c 


Fe I 


0.02 


0.22 


-0.12 


0.11 


-0.20 


0.10 


-0.25 


0.18 


0.01 


0.11 


0.00 


0.13 


Fe II 


0.07 


0.16 


-0.10 


0.10 


-0.20 


0.13 


-0.27 


0.21 


0.06 


0.10 


0.00 


0.12 


[X/Fe] 


























C 


0.58 


0.15 


0.45 


0.15 


0.46 


0.15 


0.79 


0.15 


0.44 


0.15 


0.45 


0.15 


o 


0.13 


0.15 


0.15 


0.15 


0.17 


0.15 


0.30 


0.15 


0.14 


0.15 


0.10 


0.15 


Y I 


0.78 


0.20 


0.77 


0.07 


0.85 


0.07 


1.05 


0.14 


0.49 


0.40 


0.80 


0.10 


Y II 


0.66 


0.06 


0.83 


0.05 


0.82 


0.09 


1.12 


0.09 


0.34 


0.14 


0.76 


0.20 


Zr I 


0.56 


0.06 


0.84 


0.10 


0.82 


0.08 


1.20 


0.09 


0.53 


0.11 


0.98 


0.17 


Zr II 


0.43 


0.10 


0.75 


0.07 


0.85 


0.07 


1.12 


0.07 


0.44 


0.12 


0.95 


0.07 


La II 


0.43 


0.20 


0.77 


0.06 


0.84 


0.05 


0.75 


0.09 


0.24 


0.15 


0.62 


0.09 


N(l II 


0.40 


0.10 


0.78 


O.O.'-i 


0.87 


0.05 


0.85 


0.05 


0.54 


0.06 


0.71 


0.08 


Ku 11 


0.02 


0.10 


0.00 


0.10 


0.09 


0.10 


0.0() 


0.10 


-0.07 


0.10 


0.10 


0.10 


C/O 


1.20 




0.85 




0.83 




1.32 




0.85 




0.95 




Is 


0.55 




0.79 




0.83 




1.12 




0.39 




0.87 




hs 


0.42 




0.77 




0.86 




0.80 




0.39 




0.67 




hs/ls 


-0.13 




-0.02 




0.02 




-0.32 




0.00 




-0.21 





Table 7. Dependence of derived abundances on atmospheric parameters 



Species 


(+250K) 


A log g 
(+0.3) 


A? 

(-1-0.25 km s-i) 


O I 


-0.60 


0.35 


-0.02 


Y I 


0.31 


0.02 


-0.16 


Y II 


-0.09 


0.15 


-0.16 


Zr I 


0.32 


0.03 


-0.31 


Zr II 


-0.13 


0.15 


-0.11 


La II 


0.02 


0.15 


-0.24 


Nd II 


0.01 


0.14 


-0.18 


Eu II 


-0.05 


0.15 


-0.07 
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Fig. 1. — Relationship between [Fe/H] and effective temperature for field stars. 
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Fig. 2. — Distribution of C/0 ratio for the field stars studied. 
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Fig. 3. — Examples of line synthesis used to ascertain C and O abundances, a) HD49331 
with C abundances of [C/Fe]=0.00 (short dashed), 0.05 (dotted), 0.15 (long dashed) and 0.30 
(solid) and b) HD64332 with C abundances of [C/Fe]=0.00, 0.50, 0.65 and 0.80; c) HD49331 
and d) HD64332 with O abundances of [O/Fe]=0.00, 0.10, 0.25 and 0.40. 
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Fig. 4. — Examples of line synthesis used to ascertain Zr and Nd abundances, a) HD49331 
and b) HD64332 with Zr abundances of no Zr (short dashed), [Zr/Fe]=0.00 (dotted), 0.50 
(long dashed) and 1.00 (sohd); c) HD49331 and d) HD64332 with Nd abundances of no Nd, 
[Nd/Fe]=0.00, 0.50 and 1.00. 
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6125 6130 6135 6140 6145 



Wavelength (A) 

Fig. 5.— Spectra of a K giant (top) (HD27371, Te//=5200K, log g=3.0, [Fe/H]=0.10), and 
an M (HD49331), S (HD49368) and C (HD30959) star. The increase in the strength of the 
Zr hnes is obvious. 
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Fig. 6. — Spectra of the same stars shown in Figure [5] (top to bottom: HD27371, HD49331, 
HD49368 and HD30959). The increase in the strength of the Nd hnes is obvious. 
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Fig. 7. — Trends in [Is/Fe], [hs/Fe] and [hs/ls] for different spectral types. 



-29- 



I 

9 



o 



o 



CD 



' — ID 



CD 



o 



• [Is/Fe] 



° [hs/Fe] 







0.5 

[C/Fe] 



Fig. 8. — Trends in [Is/Fe], [hs/Fe] as the [C/Fe] value increases. 
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Fig. 9. — Trends in [Is/Fe], [hs/Fe] and [hs/ls] as the C/0 ratio increases. Results from 
previous studies are also shown (squares: Smith and Lambert, 1985; triangles: Abia et al., 
2002). 



